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(54) Estimation of frequency offset, for use with short data bursts 



(57) Frequency offset estimation, suitable for use 
with bursts of up to 200 symbols, is performed in two 
stages. Firstly, an initial estimate is made, using the pre- 



amble. Secondly, during the remainder of the burst, 
those decisions which are deemed unreliable are not 
used in the estimation. 
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Description 

Background 

[0001] Figs. 1 A and 1 B illustrate a prior-art Hybrid Fib- 
er-Coax (HFC) cable system 100 that is compatible with 
the cable industry standard Data over Cable System In- 
terface Specification (DOCSIS) for providing Internet 
access to cable customers via so called Cable Modems 
(CMs). Fig. 1 A is a top-level view of the cable system. 
Fig. 1B provides additional detail of the Customer 
Premises Equipment (CPE) of Fig. 1A. In Fig. 1B, CM 
4000 provides a computer industry standard Ethernet 
interface to PC 5000 and bridges the Ethernet interface 
with the coax distribution of the cable system. CM 4000 
internally implements both an RF Modulator and an RF 
Demodulator for communications over the coax in ac- 
cordance with the DOCSIS standard. 
[0002] An RF Modulator 3000 and RF Demodulator 
1000, complementary to those of the cable modem, are 
implemented in a DOCSIS compatible Cable Modem 
Termination System (CMTS) 500, which as the name 
implies, provides termination for the Cable Modem of 
the CPE. Multiple instances of Modulator 3000 and De- 
modulator 1000 are provisioned to support all custom- 
ers with CM service. Control, MAC, Framing 2000 bridg- 
es all of the provisioned DOCSIS RF interfaces with one 
or more packet-based networks. These packet net- 
works may include local area networks, intranets, and 
the Internet. While Fig. 1 A shows the CMTS 500 imple- 
mented in a Head End or Primary Hub, theoretically it is 
possible to implement the CMTS anywhere upstream 
from the CM. Each demodulator 1000 provides outputs 
to the Control, MAC, Framing 2000 that include Detect- 
ed Symbols 1200 and a Frequency Offset Estimate 
1300. 

[0003] While the transmitter of the CM upstream mod- 
ulator and the complementary receiver in the CMTS de- 
modulator are theoretically provisioned for identical fre- 
quencies, frequency offsets may occur for a variety of 
reasons. These reasons include errors in the local os- 
cillators of the CM modulator (used to transmit the signal 
burst at the provisioned frequency) or CMTS front-end 
(used to down-convert the received signal burst to base- 
band) as well as the use of up- and down-conversion 
processes applied to the upstream channels over the 
path from CM to CMTS. These conversion processes 
are used to combine and split multiple channels over 
common communication paths for economic transport. 
Depending on severity and other system issues, fre- 
quency offsets could degrade Bit Error Rates (BERs), 
increase latencies, or completely prevent signal cap- 
ture. However, these problems may be avoided by using 
closed loop methods to minimize or eliminate the off- 
sets. 

[0004] In CMTS applications, system operation may 
be conceptually divided into normal data traffic condi- 
tions (traffic mode) and so-called ranging periods (rang- 



ing mode). Ranging is a closed loop process by which 
the CMTS manages the timing, power level, offset fre- 
quency, and equalization for the transmitter of each CM. 
Ranging is performed whenever a CM is initialized and 
5 registered by the network and also periodically (at reg- 
ular time intervals) to update its calibration. The ranging 
calibration process is performed for every CM on the 
channel and enables the system to smoothly operate at 
high effective throughput during traffic mode. 
10 [0005] As applicable to the cable system of Fig. 1 A, 
during ranging, a Frequency Offset Estimate 1300 is 
made in Demodulator 1000 of the CMTS and provided 
to the Control, MAC, and Framing 2000. The control 
functions of this block then send one or more commands 
15 downstream to the corresponding CM, remotely adjust- 
ing its transmitter frequency until the frequency offset 
observed at Demodulator 1 000 is reduced to within pre- 
defined acceptable limits. Once set during ranging, the 
offset frequency adjust continues to be used by the CM 
20 during subsequent traffic mode operation. The ranging 
determined offset frequency adjust thus effectively elim- 
inates any frequency offset observed by the CMTS, re- 
gardless of the source or sources of the offset. 
[0006] Ranging represents the most problematic op- 
25 erating condition for determining the frequency offset 
estimate, as the frequency offset is generally large dur- 
ing ranging and other CM characteristics are not yet cal- 
ibrated for optimal performance. Accordingly, symbol er- 
rors may be frequent. During traffic mode, the CM is op- 
30 erating with a much smaller effective frequency offset 
due to the compensating offset frequency adjust as- 
signed by the CMTS during ranging and overall CM op- 
eration is optimally calibrated. Symbol errors are signif- 
icantly reduced. 

35 

Demodulator Operation 

[0007] Fig. 2A provides a general conceptual block di- 
agram of the digital burst Demodulator 1000 in the 
40 CMTS 500. Front-End 600 isolates one modulated car- 
rier from the carrier multiplex in the Received Spectrum 
1100, baseband converts the signal, and passes the re- 
sulting signal 1105 to the Burst and Timing Synchroni- 
zation circuit 1500. (In other contexts the Front-End 600 
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of, the demodulator.) The Recovered Signal Samples 
1106, at the output of circuit 1500, are discrete signal 
samples at the symbol rate. 

[0008] Equalizer 1 600 compensates for signal distor- 
so tion not compensated by a pre-equalizer in the cable 
modem (CM) and also suppresses ingress noise. At the 
output of this stage, the Equalized Signal Samples r k 
(1107) are given by: 



55 /<p* 

r k =a k e + w k , 



where a^is the transmitted symbol and w^is the additive 
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noise. These samples are not yet synchronized in terms 
of carrier phase. 

[0009] Carrier Phase Synchronized Samples 1108 
are produced by the Rotator 1700 in conjunction with 
Phase Estimator 1 900. The carrier phase y k is given by: 

9^= 2knAfT s + <p 0 . 

where Af is the frequency offset and T s is the symbol 
period. 

[0010] Let $/, be the estimated carrier phase used for 
symbol k. The equalized signal r k is multiplied by 



known symbols) and a short data block consisting of 
symbols extracted by the local detector. Due to noise 
and other factors, generally including the existence of 
frequency offsets, the output of the detector will not be 

5 free of symbol errors. 

[0014] In view of the above, it is useful to make a de- 
tailed examination of the performance of the prior art fre- 
quency estimator based on the autocorrelation function 
of the signal. As shown in Fig. 2B, the prior art autocor- 

10 relation estimator bases the frequency estimation on the 
computation of 8* = PfP * By proper normalization of 
the time average of 8^ the frequency can be estimated 
as: 



15 



within Rotator 1700 so as to decide symbol a k by De- 20 
tector 1800. The decision on symbol a k is denoted d k , 
which is output as Detected Symbols 1200. 
[001 1] The quantity p^r k d is subsequently comput- 
ed. For CMTS and other applications having a high val- 
ue of the frequency drift, quasi-differential demodulation 25 
is preferred, and the carrier phase recovery algorithm 
only uses the last value of p k to determine the carrier 
phase. In such cases, the estimated carrier phase can 
thus be written as: 

30 
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Prior Art Frequency Estimators 



35 



[0012] Frequency Offset Estimator 8000 provides 
Frequency Offset Estimate 1300, an estimate of the ob- 
served frequency offset of the received signal. Prior art 
frequency estimation techniques for general burst de- 
modulator applications have used phase slope compu- 40 
tations and related time averaging, or have been based 
on the autocorrelation function of the signal. Fig. 2B il- 
lustrates a particular prior art Frequency Offset Estima- 
tor (8000- PA- AC) based on autocorrelation of the signal. 

Ci« 00 Hji i^trot^e q r»orti^i ilar nrinr art Frennennv Offset 45 

Estimator (8000-PA-PS) based on phase slope compu- 
tations. Fig. 2D provides additional detail of the Phase 
Slope Computation 8100, of Fig. 2C. 
[0013] The prior art frequency estimation techniques 
provide satisfactory precision for more general burst de- 50 
modulator applications that use long bursts of known 
symbols (such astound in long message preambles) or 
admit to time averaging over multiple consecutive 
bursts. Unfortunately, both isolated short bursts and 
symbol errors are characteristics of CMTS applications, 55 
particularly during ranging. CMTS applications employ 
isolated short ranging bursts of approximately 200 sym- 
bols, consisting of a very short preamble (typically 20-30 



N 



V*=i+i j 
2nT s L 



where N designates the number of symbols in the 
burst at hand, and L is a delay parameter defined in the 
Correlator 8620 in Fig. 2B. 

[0015] Suppose that a symbol error occurs at the k 0 
decision instant. With the assumed QPSK signal format, 
this leads to an instantaneous phase error of ±n/2. Fur- 
thermore, as quasi-differential demodulation is used, 
this ±ti/2 phase-error appears in all subsequent symbol 
decisions, and also in all p k values computed after k 0 . 
[0016] Thus, the phase of p k for the theoretical fre- 
quency estimator has a constant slope except for one 
discontinuity of ±nJ2 when the error occurs. Conse- 
quently, the phase of b k (recall that 8^ = Ptf> k . ) is cen- 
tered around the desired value, except for Lvalues with 
an error of ±nl2. The average phase of § k over the burst 
size is then shifted by ±L(n/2)/(N-L). Since the normal- 
ized frequency estimate ArT s is obtained by dividing this 
average phase by 2nL, it clearly contains an error of 
±0.25/(AM_). 

[0017] Consider a particular quantitative example us- 
ing the foregoing analysis. Suppose that we want to 
have an estimation precision of 10" 4 for a burst length 
of 200. Should a single decision error occur during this 
burst, this error would result in a normalized frequency 
estimation error of approximately 0.25/200, or approxi- 
mately 10" 3 . This is an order of magnitude greater than 
the desired accuracy of 1 0' 4 . Those skilled in the art will 
see that a similar analysis, with similar results, can be 
made for a phase slope based estimator. 
[0018] The prior art frequency estimation approaches 
are clearly not well suited to CMTS applications, as it 
has been seen that even a single symbol error results 
in substantial errors. Improved frequency estimation 
methods and devices are needed that provide precise 
frequency estimates for short bursts in the presence of 
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symbol errors. New frequency estimation approaches 
are needed that can be performed on one single burst, 
without the need for averaging the results of several 
consecutive bursts. In particular, improved frequency 
estimation techniques are needed that offer superior 
performance for CMTS applications. 

Brief Description of Drawings 

[0019] 

Figs. 1 A and 1 B illustrate a prior-art HFC cable sys- 
tem. Fig. 1 A is a top-level view of the cable system. 
Fig. 1 B provides additional detail of the CPE of Fig. 
1A. 

Fig. 2A provides internal architectural detail of De- 
modulator 1000 of Fig. 1A. 
Fig. 2B provides detail of a prior art Frequency Off- 
set Estimator 8000-PA-AC, based on Autocorrela- 
tion of the signal, and having application in the De- 
modulator 1000 of Fig. 2A. 
Fig. 2C provides detail of a prior art Frequency Off- 
set Estimator 8000-PA-PS, based on Phase Slope 
Computation, and having application in the Demod- 
ulator 1000 of Fig. 2A. 

Fig. 2D provides additional detail of a prior art 
Phase Slope Computation Function 8100, as used 
in Fig 2C. 

Fig. 3 provides detail of a first illustrative embodi- 
ment (8000-AC) of a Frequency Offset Estimator, 
based on Autocorrelation, in accordance with the 
present invention, and having application in the De- 
modulator of Fig. 2A. 

Fig. 4 provides detail of a second illustrative embod- 
iment (8000-PS) of a Frequency Offset Estimator, 
based on Phase Slope Computation, in accordance 
with the present invention, and having application 
in the Demodulator of Fig. 2A. 
Fig. 5 provides detail of an illustrative embodiment 
of an Atypical Samples Filter 9000, suitable for ap- 
plication as either the Discard Atypical Correlation 
Samples Function 8630 of Fig. 3, or the Discard 
Atypical Phase Slope Samples Function 8200 of 
Fig. 4. 

Summary 

[0020] The present invention provides methods and 
apparatus for generating precise frequency estimates 
for applications with both short bursts and symbol er- 
rors. The invention is practical, easy to implement, ro- 
bust against decision errors, and requires only one burst 
to achieve frequency estimation with high precision 
(around 10* 4 in AfT s in the illustrative embodiments). 
[0021] The invention performs frequency estimation 
over both the burst preamble, during which known sym- 
bols are transmitted, and also during the burst's data 
packet, which is subsequent to the preamble and ex- 



tracted by the local detector. During the preamble, an 
initial frequency estimate is obtained. This estimate is 
based on a time average of either phase or correlation 
samples. Atypical phase or correlation samples, attrib- 
5 utable to detector symbol errors during the data packet, 
are detected and filtered, so as to avoid including the 
atypical samples in a time-average used to provide the 
frequency estimate. 

[0022] In a first embodiment correlation samples are 
10 time averaged, and atypical correlation samples are 
suppressed prior to correlation time averaging. In a sec- 
ond embodiment, phase slope values are time aver- 
aged, and atypical values of phase slope are sup- 
pressed prior to phase slope time averaging. 

15 

Detailed Description 
Atypical Sample Suppression 

20 [0023] Recall that in the presence of a single symbol 
error, the phase of p k for the prior art autocorrelation fre- 
quency estimator of Fig. 2B has a constant slope except 
for one discontinuity of ±rc/2 associated with the symbol 
error. Consequently, the phase of 5^ is centered around 

25 the desired value, except for L values with an error of 
±kJ2. The present invention, illustrated in Figs. 3-5, im- 
proves upon the prior art estimator by eliminating the 
secondary peak via appropriate non-linear filtering. 
More specifically, the aim of this filtering is to eliminate 

30 from the integration on 5* the values known to be wrong. 
These values are those whose argument differs from 
that of the current average of 8^ by more than a certain 
angle v. 

[0024] In preferred embodiments a preamble of 26 

35 known symbols is used. Thanks to the preamble, it is 
possible to have a first coarse frequency estimate at the 
beginning of the data symbols. Using this value, it is pos- 
sible to determine for each new b k whether it should be 
kept or not in the sum. If it is, this new value of the sum 

40 is used to filter the next symbols. Thus, if the coarse 
estimate provided by the preamble is precise enough, 
the secondary peaks are totally eradicated. Simulations 
of the illustrative embodiments have confirmed that the 
precision of the coarse frequency estimate provided by 

45 the preamble does support elimination of the secondary 
peaks. Those skilled in the art will appreciate that while 
the above description has focused on the elimination of 
atypical correlation samples for an autocorrelation fre- 
quency estimator, the invention is equally applicable to 

so elimination of atypical phase slope samples for a phase 
slope computation frequency estimator. 

Autocorrelation Embodiment 

55 [0025] Fig. 3 shows an Autocorrelation Frequency Es- 
timator (8000-AC) in accordance with the present inven- 
tion. As the name suggests, the frequency estimation 
algorithm is based on the computation of Correlation 



4 



EP 1 235 401 A1 



8 



Samples S k = p^ L (8625), which ar generated by 
Correlator 8620. 

[0026] The discard atypical values and time averag- 
ing functions (8630 and 8640, respectively) operate on 
an exponential function involving the phase slope. At the 
output of the Correlation Time Averaging stage 8640, 
we have p k expG27tLAf k T) 8645, where atypical values 
have been suppressed by Discard Atypical Correlation 
Samples 8630. Discard Atypical Correlation Samples 
8630 acts such that 5^ has no effect on the running time 
average whenever the absolute value of arg[8^ exp(- 
j2izLAf k T)] is greater than a fixed threshold v. 
[0027] Arg(.) 8650 delivers the phase 27tLAf k T, which 
is used as unnormalized frequency offset estimate 
1300-U. Normalized frequency offset estimate 1300-N 
may optionally be provided as output. 

Phase Slope Computation Embodiment 



Correlation Samples Function 8630 of Fig. 3, or the Dis- 
card Atypical Phase Slope Samples Function 8200 of 
Fig. 4. This filter acts such that the current sample 901 0 
has no effect on the running time average of samples 

5 whenever the absolute value of the difference between 
the time average of samples 9020 and current sample 
9010 is greater than the typical sample variance 9040. 
Block 9200 compares Current Sample 9010 with Time 
Average Of Samples 9020 to generate Current Sample 

10 Variance 9250. This in turn is compared against prede- 
fined Typical Sample Variance 9040 to produce Current 
Sample is Typical 9350, a logical binary valued control 
output. This control can be used directly by a subse- 
quent time averaging circuit to qualify/suppress the use 

15 of the current sample. Alternatively, or in combination 
with direct use, this control can be used to zero the Filter 
Sample value 9060 provided to the time averaging cir- 
cuit. 



[0028] Figs. 4 shows a Phase Slope Computation 20 
Frequency Estimator (8000-PS) in accordance with the 
present invention. The Phase Slope Computation 8100 
of Fig. 2D, is preferably used as the first stage of the 
improved Frequency Estimator (8000-PS) of Fig. 4. 
Within said Phase Slope Computation 8100, the re- 25 
ceived signal r k 1 1 07 is multiplied (via 81 05) by the com- 
plex conjugate of the detected symbol d k 1200. The 
function Arg(.) (i.e., the phase calculation function) 81 20 
gives the phase difference 8110 between the received 
and detected signals. Slope Computation 8140 com- 30 
putes the phase slope 81 50 from two phase-values sep- 
arated by L samples. 

[0029] Phase Slope Time Averaging 8300 outputs an 
estimate of the frequency offset 2izAf k T (1300-U) using 
the phase slope values computed up to time k. This un- 35 
normalized frequency offset, or the normalized frequen- 
cy offset estimate Af k (1 300-N), is updated at the symbol 
rate. 

[0030] Discard Atypical Phase Slope Samples 8200 
acts to discard the instantaneous phase slope values 40 
8150 which differ from the current unnormalized fre- 
quency offset estimate 2tiA^T(1 300-U) by more than a 
fixed threshold u.. If the difference between the two sig- 
nals exceeds u., then this circuit decides that a symbol 
error has occurred and inai the current slope value 45 
should not be included in the time-averaging function. 
That is, X k has no effect on the running time average 
whenever the absolute value of 2nAf k T-X k is greater 
thanu,. Eitherthe unnormalized frequency offset 1300-U 
or normalized frequency offset estimate 1 300-N may be so 
provided as output to the Control, MAC, Framing 2000 
block of CMTS 500. 



Atypical Samples Filter 

[0031] Fig. 5 provides detail of a conceptual illustra- 
tive embodiment of an Atypical Samples Filter 9000, 
suitable for application as either the Discard Atypical 



55 



Implementation Considerations 

[0032] In comparing the two illustrative embodiments, 
the phase function on which the phase slope based es- 
timator of Fig. 4 operates is a "wrapped" function that is 
periodic between -n and +71. The significance of this is 
that boundary tests and conditional calculations must be 
performed to "unwrap" this function prior to the compu- 
tation of the phase slope and thereby the frequency off- 
set. The correlation-based estimator of Fig. 3 has no 
need to perform the extra phase unwrapping operations, 
and thus may offer implementation advantages. 
[0033] In considering the implementation boundary 
between hardware (H/W) and software (S/W), remem- 
ber that the Demodulator 1000 of Fig. 1A does not use 
the frequency estimation result. Instead it is sent by the 
CMTS to the CM as an offset-frequency-adjust to control 
the CM transm it-frequency. This is not a real time pro- 
cedure, and only one offset frequency adjust is to be 
sent per burst. However, the frequency error estimate is 
extracted from a real-time sampled signal with one sam- 
ple per symbol. In our preferred implementation, all 
functions are implemented in hardware except the 
phase calculation in box 8650 and optional normaliza- 
tion in box 8660 of Fig. 3, and the optional normalization 
in h 0 v S400 in Fig. 4, The reason is that those compu- 
tations are performed once per burst, which means that 
they do not need high-speed electronics. 
[0034] A fully software implementation would consist 
of sending the equalized samples r k and detected sym- 
bols d k to the software and make all computations in 
software. While such an implementation is within the 
scope of the present invention, it is not likely to be an 
efficient implementation. In general, real-time computa- 
tions (at the symbol rate) are best made in hardware 
and low-speed computations (at the burst rate) are best 
made in software. 
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Conclusion 

[0035] Although the present invention has been de- 
scribed using particular illustrative embodiments, it will 
be understood that many variations in construction, ar- 
rangement and use are possible consistent with the 
teachings and within the scope of the invention. For ex- 
ample, bit-widths, clock speeds, and the type of tech- 
nology used may generally be varied in each component 
block of the invention. Also, unless specifically stated to 
the contrary, the value ranges specified, the maximum 
and minimum values used, or other particular specifica- 
tions (such as the desired accuracy of the frequency es- 
timate), are merely those of the illustrative or preferred 
embodiments, can be expected to track improvements 
and changes in implementation technology, and should 
not be construed as limitations of the invention. Func- 
tionally equivalent techniques known to those skilled in 
the art may be employed instead of those illustrated to 
implement various components or sub-systems. It is al- 
so understood that many design functional aspects may 
be carried out in either hardware (i.e., generally dedi- 
cated circuitry) or software (i.e., via some manner of pro- 
grammed controller or processor), as a function of im- 
plementation dependent design constraints and the 
technology trends of faster processing (which facilitates 
migration of functions previously in hardware into soft- 
ware) and higher integration density (which facilitates 
migration of functions previously in software into hard- 
ware). 

[0036] Specific variations within the scope of the in- 
vention include, but are not limited to: the use of any of 
a variety of techniques for identifying atypical samples 
and for suppressing their impact on sample averages, 
the use of any of a variety of preamble lengths and not 
just via the preferred length of 26 symbols, and the use 
of either or both of the unnormalized and normalized fre- 
quency offset. 

[0037] All such variations in design comprise insub- 
stantial changes over the teachings conveyed by the il- 
lustrative embodiments. The names given to intercon- 
nect and logic are illustrative, and should not be con- 
strued as limiting the invention. It is also understood that 
the invention has broad applicability to other communi- 
cations and network applications, and is not limited to 
the particular application or industry of the illustrated 
embodiments. The present invention is thus to be con- 
strued as including all possible modifications and varia- 
tions encompassed within the scope of the appended 
claims. 
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Claims 

1 . A subsystem for estimating the frequency offset be- 
tween a transmitter and a receiver, the receiver hav- 
ing signals including pre-detection signal samples 
and detected symbols, the subsystem comprising: 



a) a first processing stage, coupled to the pre- 
detection signal samples and the detected 
symbols, the first processing stage generating 
first processed samples; 

b) a time averaging stage, having a samples in- 
put and a first estimate indication output, the 
first estimate indication being a time average of 
the input samples; and 

c) an atypical samples stage, coupled to re- 
ceive the first processed samples and the first 
estimate indication, the atypical samples stage 
having a sample status indication selectively 
characterizing particular first processed sam- 
ples as typical. 

2. The subsystem of claim 1 , further wherein the sam- 
ple status indication is coupled to the time averag- 
ing stage to suppress inclusion of atypical samples. 

20 3. The subsystem of claim 2, further wherein the sam- 
ples input of the time average stage is coupled to 
an output of the atypical samples stage. 

4. The subsystem of claim 1 , further wherein the atyp- 
25 ical samples stage provides filtered samples to the 
samples input of the time average stage, the sam- 
ples being generated from the first processed sam- 
ples by suppressing atypical samples. 

30 5. The subsystem of claim 4, wherein zero samples 
are substituted for the atypical samples to perform 
the suppression. 

6. The subsystem of claim 5, wherein the substitution 
35 is performed using a signal multiplexor. 

7. The subsystem of claim 6, wherein the multiplexor 
is controlled by the signal status indication. 

40 8. The subsystem of claim 1 , wherein the atypical 
samples stage creates an intermediate current 
sample variance. 

9. The subsystem of claim 1 , wherein the atypical 
45 samples stage performs first and second compari- 
son operatioTisrwherein the first comparison oper- 
ation evaluates the current sample against the av- 
erage of previous samples and creates a current 
sample variance, and the second comparison op- 

50 eration evaluates the current sample variance 
against a typical sample variance. 

1 0. The subsystem of claim 1 , wherein the first process- 
ing stage includes a correlator and the first proc- 

55 essed samples are correlation samples. 



11. The subsystem of claim 10, further including a 
phase calculation stage having an input coupled to 
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the first estimate indication and generating a sec- 
ond estimate indication, the second estimate indi- 
cation providing an unnormalized frequency offset 
estimate. 

12. The subsystem of claim 11 , further including a nor- 
malization stage having an input coupled to the sec- 
ond estimate indication and generating a third esti- 
mate indication, the third estimate indication provid- 
ing a normalized frequency offset estimate. 

13 The subsystem of claim 1, wherein the first process- 
ing stage performs a phase slope computation, the 
first processed samples are phase slope samples, 
and the first estimate indication provides an unnor- 
malized frequency offset estimate. 

1 4 The subsystem of claim 1 3, further including a nor- 
malization stage having an input coupled to the first 
estimate indication and generating a second esti- 
mate indication, the second estimate indication pro- 
viding a normalized frequency offset estimate. 

15 The subsystem of claim 13, wherein the phase 
slope computation includes a phase calculation 
sub-stage and a slope computation sub-stage. 

16 The subsystem of claim 1 , wherein the subsystem 
is operated over both the burst preamble and during 
the burst's data packet, and wherein the atypical 
samples stage does not suppress samples during 
at least part of the burst preamble. 

17. The subsystem of claim 16, wherein the preamble 
uses 26 known symbols. 

18. The subsystem of claim 16, wherein the atypical 
samples stage does not suppress samples over the 
entire burst preamble. 

1 9 A method of estimating a frequency offset between 
a transmitter and a receiver based on the received 
spectrum at the receiver, the transmitter sending 
the receiver bursts having a preamble and data, the 
method comprising: 

a) providing processed samples derived from 
the received spectrum, including a current sam- 
ple; 

b) providing a typical sample variance; 

c) generating a time average of the processed 
samples; and 

d) during the burst data, excluding the current 
sample in the time average when the difference 
between the current sample and the time aver- 
age exceeds the typical sample variance. 

20. The method of claim 19, further including: 



during the burst preamble, including the current 
sample in the time average. 

21 . The method of claim 20, wherein the preamble uses 
5 26 known symbols. 

22. The method of claim 20, wherein the processed 
samples are correlation samples. 

w 23. The method of claim 22, further including calculat- 
ing the phase of the time average to provide an un- 
normalized frequency offset estimate. 

24 The method of claim 23, further including perform- 
ing normalization to provide a normalized frequency 
offset estimate. 



25. The method of claim 20, wherein the processes 
samples are phase slope samples and the time av- 
erage provides an unnormalized frequency offset 
estimate. 



26 The method of claim 25, further including perform- 
ing normalization to provide a normalized frequency 

25 offset estimate. 

27 A burst demodulator for a receiver, the receiver re- 
ceiving bursts from a transmitter, the burst demod- 
ulator comprising: 

a) a stage generating pre-detection signal sam- 
ples; 

b) a detector stage generating detected sym- 
bols; and 

c) a subsystem for estimating the frequency off- 
set between the transmitter and the receiver, 
the subsystem including 

i. a first processing stage, coupled to the 
pre-detection signal samples and the de- 
tected symbols, the first processing stage 
generating first processed samples; 

ii. a time averaging stage, having a sam- 
ples input and a first estimate indication 

. . «^timQta inHiratinn heina a 

OUipUL, Liie MiSi colmiimiw ^ 

time average of the input samples; and 
Hi. an atypical samples stage, coupled to 
receive the first processed samples and 
the first estimate indication, the atypical 
samples stage having a sample status in- 
dication selectively characterizing particu- 
lar first processed samples as typical. 

28. The burst demodulator of claim 27, wherein the first 
processing stage includes a correlator and the first 
processed samples are correlation samples. 

29. The burst demodulator of claim 28, further including 
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a phase calculation stage having an input coupled 
to the first estimate indication and generating a sec- 
ond estimate indication, the second estimate indi- 
cation providing an unnormalized frequency offset 
estimate. 

30 The burst demodulator of claim 29, further including 
a normalization stage having an input coupled to the 
second estimate indication and generating a third 
estimate indication, the third estimate indication 
providing a normalized frequency offset estimate. 

31 The burst demodulator of claim 27, wherein the first 
processing stage performs a phase slope compu- 
tation the first processed samples are phase slope 
samples, and the first estimate indication provides 
an unnormalized frequency offset estimate. 

32 The burst demodulator of claim 31 , further including 
a normalization stage having an input coupled to the 
first estimate indication and generating a second 
estimate indication, the second estimate indication 
providing a normalized frequency offset estimate. 

33 The burst demodulator of claim 31, wherein the 
' phase slope computation includes a phase calcula- 
tion sub-stage and a slope computation sub-stage. 

34 The burst demodulator of claim 27, wherein the sub- 
system is operated over both the burst preamble 
and during the burst's data packet, and wherein the 
atypical samples stage does not suppress samples 
during at least part of the burst preamble. 

35. The burst demodulator of claim 34, wherein the pre- 
amble uses 26 known symbols. 

36 The burst demodulator of claim 34, wherein the 
atypical samples stage does not suppress samples 
over the entire burst preamble. 

37 A cable modem termination system (CMTS) for a 
cable system, the cable system including at least 
one cable modem having a transmitter, the CMTS 
comprising: 
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i. a first processing stage, coupled to the 
pre-detection signal samples and the de- 
tected symbols, the first processing stage 
generating first processed samples; 

ii. a time averaging stage, having a sam- 
ples input and a first estimate indication 
output, the first estimate indication being a 
time average of the input samples; and 
Hi. an atypical samples stage, coupled to 
receive the first processed samples and 
the first estimate indication, the atypical 
samples stage having a sample status in- 
dication selectively characterizing particu- 
lar first processed samples as typical. 

38 The CMTS of claim 37, wherein the first processing 
stage includes a correlator and the first processed 
samples are correlation samples. 

20 39 The CMTS of claim 38, further including a phase 
calculation stage having an input coupled to the first 
estimate indication and generating a second esti- 
mate indication, the second estimate indication pro- 
viding an unnormalized frequency offset estimate. 

25 40 The CMTS of claim 39, further including a normali- 
zation stage having an input coupled to the second 
estimate indication and generating a third estimate 
indication, the third estimate indication providing a 
30 normalized frequency offset estimate. 

41 The CMTS of claim 37, wherein the first processing 
stage performs a phase slope computation, the first 
processed samples are phase slope samples, and 

35 the first estimate indication provides an unnormal- 
ized frequency offset estimate. 

42 The CMTS of claim 41 , further including a normali- 
zation stage having an input coupled to the first es- 

40 timate indication and generating a second estimate 
indication, the second estimate indication providing 
a normalized frequency offset estimate. 
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a) a modulator; 

b) a network interface including control, media- 
access-control, and framing functions; and 

c) ademodulatorcoupledto receive bursts from 
the cable modem transmitter, the demodulator 
including a stage generating pre-detection sig- 
nal samples, a detector stage generating de- 
tected symbols, and a subsystem for estimating 
the frequency offset between the cable modem 
transmitter and the demodulator, the subsys- 
tem including 
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43. The CMTS of claim 41, wherein the phase slope 

computation mciuues a H" a<> ^ w ~ 

stage and a slope computation sub-stage. 

44 The CMTS of claim 27, wherein the subsystem is 
operated over both the burst preamble and during 
the burst's data packet, and wherein the atypical 
samples stage does not suppress samples during 
at least part of the burst preamble. 

45. The CMTS of claim 44, wherein the preamble uses 
26 known symbols. 

46. The CMTS of claim 44, wherein the atypical sam- 
ples stage does not suppress samples over the en- 
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tire burst preamble. 



47. A method of adjusting the frequency of a cable mo- 
dem, the cable modem having a transmitter that 
sends bursts to the demodulator of a cable modem 5 
termination system, the bursts having a preamble 
and data, the demodulator having a received spec- 
trum at its input, the method comprising: 

a) commanding the cable modem to operate in « 
ranging mode; 

b) providing processed samples derived from 
the received spectrum, the particular proc- 
essed sample at the current time being referred 
to as the current sample; 

c) defining a typical sample variance and an ac- 
ceptable frequency offset limit; 

d) generating a running time average of the 
processed samples; 

e) during the burst preamble, including the cur- 
rent sample in the time average; 

f) during the burst data, excluding the current 
sample in the time average when the difference 
between the current sample and the time aver- 
age exceeds the typical sample variance; 

g) deriving a frequency offset estimate from the 
time average; and 

h) commanding the cable modem to adjust its 
transmit frequency until thefrequency offset es- 
timate is reduced to within the frequency offset 
limit. 

48. The method of claim 47, wherein the preamble uses 
26 known symbols. 

35 

49. The method of claim 47, wherein the processed 
samples are correlation samples. 

50. The method of claim 49, wherein the derivation of 
the frequency offset estimate includes calculating 
the phase of the time average. 

51. The method of claim 50, wherein the derivation of 
the frequency offset estimate further includes per- 
forming normalization on the phase calculation re- « 
suit. 

52. The method of claim 47, wherein the processes 
samples are phase slope samples. 
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53. The method of claim 52, wherein the derivation of 
the frequency offset estimate includes performing 
normalization on the time average. 
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